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We perform a perturbative calculation of the influence of dynamical HISQ fermions on the per- 
turbative improvement of the gluonic action in the same way as we have previously done for asqtad 
fermions. We find the fermionic contributions to the radiative corrections in the Liischer-Weisz 
gauge action to be somewhat larger for HISQ fermions than for asqtad. Using one-loop pertur- 
bation theory as a test, we estimate that omission of the fermion-induced radiative corrections in 
dynamical asqtad simulations will give a measurable effect. The one-loop result gives a systematic 
shift of about —0.6% in ri on the coarsest asqtad improved staggered ensembles. This is the correct 
sign and magnitude to explain the scaling violations seen in "I>s = /bV-^b on dynamical lattice 
ensembles. 
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I. INTRODUCTION 

The Fermilab, MILC, HPQCD and UKQCD Collabo- 
rations are involved in an ambitious programme of high 
precision predictions of phcnomcnologically relevant pa- 
rameters from QCD using unquenched lattice simula- 
tions pQ. 

Central to this programme is the perturbative improve- 
ment of the fermionic and gluonic action and operators to 
remove significant sources of scaling violation in the lat- 
tice simulation results. This body of work is based on the 
Symanzik-improved staggered-quark formalism, specifi- 
cally the use of the asqtad action [2J. More recently, 
the Highly Improved Staggered Quark (HISQ) action has 
been used to further suppress taste-changing interactions 
and to allow the use of heavier quarks at the same lattice 
spacing by removing tree-level C((ma) 4 ) artifacts from 
the valence quark action [3J. 

To maintain the same level of improvement when these 
actions are used to describe the sea quarks |U [S], we 
should include the effect of fermion loops on the radia- 
tive terms in the Symanzik-improved gauge action. This 
has recently been done to 0(Nfa s a 2 ) for the asqtad ac- 
tion [S] and in this paper we update that calculation to 
include dynamical HISQ fermions instead. Some prelim- 
inary results can be found in Ref. [7]. We note that the 
corrections are larger for HISQ than for asqtad. 

In the second part of this paper, Sec. |IV| we consider 
what effect the 0(Nfa s a 2 ) will have in a practical sim- 
ulation, particularly on the scale-setting parameters f\ 
and f derived from the static quark potential. 

The MILC and UKQCD Collaborations have already 
used dynamical asqtad quarks to generate a large set of 
Monte Carlo lattice ensembles, including ones with very 
light sea quarks (MILC, e.g. [T]) and ones with large num- 
bers of independent configurations (UKQCD [8]). 



The gauge action used, however, omitted the asqtad 
0(N fOc s a 2 ) radiative improvements (they were not then 
known). It has been observed that the quantity $b = 
fsV^B shows a +2% scaling violation on the dynamical 
"coarse" ensembles (lattice spacing a ~ 0.12 fm) relative 
to the "fine" (a ~ 0.09 fm) [pj. This scaling violation was 
not seen for corresponding quenched lattices [10] . 

Unless there is a subtle (and therefore unlikely) cancel- 
lation, the quenched result suggests that the (quenched) 
gluonic and valence staggered actions are not the prob- 
lem. If the asqtad action is suitable for the valence 
quarks, it seems likely it is equally suitable for the sea 
quarks. The scaling violation is therefore argued to arise 
from the omission of the 0(Nfa s a 2 ) radiative corrections 
to the gluonic action. It is certainly plausible that the 
fermionic contributions could have such an effect; they 
are, after all, large enough to reverse the sign of some of 
the radiative couplings in the action [5]. 

In these calculations fx = r x /a has been used to set 
the scale, i.e. to convert from dimcnsionlcss lattice re- 
sults to physical predictions. We therefore attempt to 
estimate, at least semi-quantitatively, whether the omis- 
sion of the asqtad (D(Nfa s a 2 ) radiative corrections would 
have a measurable effect on the static potential and, in 
particular, on the scale-setting parameters f\ and fo used 
to convert from dimensionless lattice results to physical 
predictions. We do this using one-loop perturbation the- 
ory. We treat such a result as indicative: we do not 
rule out higher loop and non-perturbative contributions, 
but argue that if one-loop perturbation theory predicts 
a measurable result then it is likely to persist when we 
include other contributions. 

Using one-loop perturbation theory, we find that in- 
cluding the fermionic radiative corrections to the gauge 
action would lead to a 0.65% decrease in f\ on the coarse 
ensembles and no change on the fine. The sign and mag- 
nitude of these shifts are robust under reasonable varia- 
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tions in fitting parameters. This would equate to a 0.65% 
increase in a on the coarse ensembles. The quantity (f>s 
scales as a -3 / 2 , so the 0(Nfa s a 2 ) corrections would lead 
to a 1% decrease in $b on the coarse ensembles and no 
effect on the fine. 

This shift is very close to what has been observed and 
we therefore suggest that the anomalous upward shift in 
$b is in large part due to the omission of the 0(NfCt s a 2 ) 
radiative corrections to the gluon action. We therefore 
predict that other observables scaling with a similar neg- 
ative power of a should exhibit similar scaling violations 
that ought to become noticeable if these observables are 
measured to similar accuracy. 

II. ON-SHELL IMPROVEMENT 

We begin by briefly reviewing how radiative improve- 
ment works at 0{a s a 2 ). 

Starting from the Symanzik tree-level-improved gauge 
action, the Coulomb self-energy from the 1-loop radiative 
corrections is (see Section 4 of Ref. [TT] and Eq. (44) of 
Ref. H2I): 

w(k) oc k 2 +a s ( ai k 2 - (3 Q k 2 ln(k 2 ) + a 2 k^ + a 3 (k 2 ) 2 ) 

(1) 

where 

k 2 =£fc 2 , kW^t (2) 

i=l i=l 

The first two terms in the brackets are absorbed into 
the scheme definition of ay (Eq. (46) of Ref. [12]) and 
do not concern us. The last two terms in the brackets 
are lattice artifacts and are 0(a s a 2 ). It is the goal of 
radiative improvement to remove these terms. 

To improve the gauge theory at 0(a s a 2 ) we introduce 
appropriate radiative counterterms into the gauge action. 
There are four such dimension-6 counterterms: three glu- 
onic operators (named by Liischer and Weisz [T3] as "pla- 
nar rectangles" , "parallelograms" and "bent rectangles" , 
with coefficients c\ , c 2 and C3 respectively) plus the static 
quark operator, 

c 4 a 2 ¥\7 ■ , (3) 

which contributes specifically to the static-quark poten- 
tial. 

The action normalisation condition 

c + 8( Cl +c 2 ) + 16c3 = l (4) 

ensures we get the correct gauge action in the continuum 
limit and fixes cq (the coefficient of the plaquette) , given 
the other coefficients. 

On-shell observables will remain unchanged under field 
redefinitions using the equations of motion. If we confine 
our attention to on-shell quantities, we can exploit this 



to set one of the Cj to zero [37]. The usual choice is to 
set C3 = 0. 

Looking at the terms in Eqn. in more detail, the 
a 2 term breaks rotational symmetry and its effect on the 
static-quark potential is given by the Fourier transform 

The leading 0(a s a 2 ) behaviour is ~ a s a 2 /r 3 with C(a 4 ) 
corrections that break rotational symmetry 

The effect of radiative improvement on the static po- 
tential is to set a 2 = 0, which therefore also restores 
rotational invariance of the static potential (at this level) 
and gives the correct Coulomb coefficient ay [12] . 

The 03 term, in contrast, already preserves rotational 
invariance: 

^ rfik 

SV 3 (r) ~ a s a 2 J ^ j^e^ + G(a 4 ). (6) 

with the leading contribution to the static-quark poten- 
tial being the 3D Kronecker (5 r> o (seen by changing vari- 
ables to Zi — e~ hiTi ). This, as Snippe points out [T2"] . 
does not affect the potential at non-zero r = |r| and will 
therefore not contribute to the scale setting parameters 
fi and f . In general, however, we do need to remove it: 
as well as the contact term there will be an effect for r > 
at higher order (i.e. at 0(a 4 )) because the denominator 
in the Symanzik tree-level Coulomb propagator will not 
exactly cancel the k 2 from the Feynman rules owing to 
differences in their definitions. 

Both C3 and C4 contribute to the a 3 term [TTJ [T3] so, 
with C3 = fixed as above, we can only remove it by 
introducing the static quark counterterm into the theory 
|15l 116] . i.e. by choosing an appropriate, and non-zero, 
value for C4. This has the effect of introducing staples 
onto temporal Wilson lines, which must be included in 
numerical simulations. Similarly, the C4 contact term will 
be important in, for instance, the T(2S — IS) mass split- 
ting. A contact term gives a contribution proportional 
to the square of the wavefunction at the origin. This is 
clearly different for the two states concerned, and will 
change the mass splitting [35]. We will not, however, 
consider the contact term in detail in this paper. 

A. The calculation 

Contact term aside, with c 3 = we thus need to de- 
termine ci and c 2 to complete the on-shell improvement. 
Given two independent quantities Q\ and Q 2 with ex- 
pansions 

Qi = Qi + w^a) 2 + d ljCj ([ia) 2 + O {{^a) A ) , (7) 

in powers of (/za), where fi is some energy scale, we obtain 
the 0(a 2 ) matching condition 

dijCj = ~Wi ■ (8) 
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FIG. 1: A plot of the fermionic contributions to the one-loop 
A meson self-energy /m against {ma) 2 . The vanishing of 



,0) 



jm in the infinite- volume limit can be seen clearly. 



Since this equation is linear, both sides can be decom- 
posed into a gluonic and a fermionic part; the gluonic 
part is known |12l 117] and is independent of the fermion 
action. 

In this paper, we focus on the fermionic contribu- 
tion to the radiative improvement of the gluon action. 
Such contributions come from quark loops, which there- 
fore cannot change the tree-level coefficients compared to 
the quenched case [17]. To compute the one-loop HISQ 
fermionic corrections to the gluon action, we will follow 
the same procedure as in the case of the asqtad action [6] , 
using as our two quantities Qi the three-gluon coupling 
and the mass of the so-called twisted A meson [T3] . 



B. Lattice perturbation theory 

We use lattice perturbation theory to calculate the ra- 
diative corrections. The (unsmeared) link variables 
are expressed in terms of the gauge field as 



Ufj,{x) = exp (gaAfj, (x + 



(9) 



which, when expanded in powers of g, leads to a per- 
turbative expansion of the lattice action, from which the 
perturbative vertex functions can be derived. 

The gauge field A^ is Lie algebra-valued, and can be 
decomposed as 



A,{x) = Y, A » ta 



(10) 



with the t a being anti-Hermitian generators of SU(AT), 
where N = 3 in the case of QCD. 

The improved Liischer-Weisz action that we study 



S = Yl i c oPo(x) + *Pi(x) + c 2 P 2 (x)} 



(11) 



with Co 



12 ' 



c (0) 



c 2 ) = 1 and, at tree level, c, 
- 0. The terms 



Co) 



Pa — ^ Ufj, v , 

P 2= (U^a + U^au + U^u + f7 CT; _ M;i/ ) , (12) 

are made up of appropriate traced, closed contours of 
gauge links. The notation here is that /i, v and a are 
summed over positive values and negative subscripts de- 
note hermitian-conjugatcd gauge links. 

The HISQ fermionic action is defined by an iterated 
smearing procedure with reunitarisation: 



° Pu(3) ° Pi 



Fat7)[£/] (13) 

where Pu(3) denotes the polar projection onto U(3) (as 
used in simulations [5], and not SU(3)), and the Fat7 
and modified asq smearings are defined in Ref . [3] . 

To handle the complicated form of the vertices and 
propagators in lattice perturbation theory, we employ a 
number of automation methods [19 } [20 } l21 | l22 l |2% 1 [24] 
that are based on the seminal work of Liischer and Weisz 
[T7] and are implemented in the HiPPy package [Ti?ll2T)] . 

The multi-level smearing of the gauge fields employed 
in the HISQ action presents particular problems when 
deriving the Feynman rules, even when employing auto- 
mated techniques. The solution to these is discussed in 
Refs. 0120]. 

Unless otherwise stated, we shall use g 2 as the pertur- 
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bative expansion parameter (rather than a s — f-), with 
expansions written in the form: 



2 (i) 

9 q 



0(g 4 



(14) 



The goal of this paper is to determine the fermionic con- 



tributions to Ci' and c\ L> , with c\ 

Since we will only consider fermionic loops, we do not 
need to concern ourselves with the gauge fixing, Haar 
measure and Fadeev-Popov ghost terms that appear in 
the gluonic portion of the perturbative Lagrangian. 

The loop integrals of continuum perturbation theory 
are replaced by finite sums over the points of the recip- 
rocal lattice in lattice perturbation theory. We carry out 
these sums exactly rather than using a stochastic estima- 
tor. 
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C. Twisted boundary conditions 

We work on a four-dimensional Euclidean lattice of 
length La in the x and y directions and lengths L z a, L t a 



FIG. 2: Plots of Oq ' against m q a (left) and of against m q a (right) with the fits shown for comparison. 



in the z and t directions, respectively, where a is the 
lattice spacing and L,L z ,L t are even integers. In the 
following, we will employ twisted boundary conditions 
[25j for the same purpose and in essentially the same way 
as in Refs. jT21 [17]. The twisted boundary conditions 
we use for gluons and quarks are applied to the (x,y) 
directions and are given by (y = x, y) 

u^{x + Lu) = n u u tl (x)n- 1 , 

^(x + Lv) = ^(x)^ 1 , (15) 

where the quark field ty sc (x) becomes a matrix in smell- 
colour space [26] by the introduction of a new SU(iV) 
quantum number "smell" in addition to the quark colour. 
In the (z,t) directions, we apply periodic boundary con- 
ditions. 

These boundary conditions lead to a change in the 
Fourier expansion of the fields: 

z 1 k 

*a(*) = NL ! L L E V**^) (16) 
2 ' P 

where the matrices are given by (up to an arbitrary 
phase, which may be chosen for convenience) 

r fe = n^nj 1 (17) 

and in the twisted (x, y) directions the momentum sums 
are now over 

NL NL . . , . 

p v = mn u , — < n v < — , v = (x,y) , (18) 

where m = ^ . The zero modes (n x = n y = mod N) 
are omitted from the sum in the case of the gluons, but 
not the quarks. The momentum sums for quark loops 
need to be divided by N to remove the redundant smell 
factor. 



We may consider the continuum limit of the twisted 
theory as a Kaluza-Klein theory in the infinite [z, t) 
plane. Denoting n = (n x , n y ), the stable particles in the 
(z,t) continuum limit of this effective theory are called 
the A mesons (n = (1, 0) or n = (0, 1)) with mass m and 
the B mesons (n = (1,1)) with mass v2m [12] . 

D. Small-mass expansions 

Even though we are ultimately interested in the radia- 
tive corrections in the chiral limit, we cannot set m q a = 
straightaway: the correct way to approach the chiral limit 
is to maintain m q /m > C as we take m q a — > and 
ma — * 0, where C is a constant determined by the re- 
quirement that a Wick rotation can be performed with- 
out encountering a pinch singularity [6] . 

We therefore adopt the following procedure to extract 
the 0(a 2 ) lattice artifacts: First, we expand some ob- 
servable quantity Q in powers of ma at fixed m q a: 

Q(ma,m q a) = a\^\m q a) + a^\m q a)(ma) 2 + 

O ( (ma) 4 , (ma) 4 log (ma) ) (19) 

where the coefficients in the expansion arc all functions 
of m q a. There is no term at O ((ma) 2 log (ma)) since 
the gluon action is improved at tree- level to 0(a 2 ) [T2"] . 
Then, we expand the coefficients a^\m q a) in powers of 
m q a. 

For a^\m q a) we have [57] 

af\m q a) = b$ log(m q a) + a$ . (20) 

Since we expect a well-defined continuum limit, 
a^\m q a) cannot contain any negative powers of m q a, 
but, depending on the quantity Q, it may contain loga- 
rithms; 6qq^ is the anomalous dimension associated with 
Q, and can be determined by a continuum calculation. 
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There can be no terms in (m q a) for n > since there is 
no counterterm in the gluon action that can compensate 
for a scaling violation of this kind. 

AQ) 



For a*} (rriqCL) we find 



AQ) 



AQ) 



(m q a) 



AQ) 



(m q a) 2 



(4$ + 6 2 Q 2 } log(m 9 a)) (m q a) 2 + O ((m q a) 4 ) . (21) 

After multiplication by (ma) 2 , the (m q a)~ 2 contribution 

gives rise to a continuum contribution to Q, and a^-^ is 
calculable in continuum perturbation theory. There can 
be no term in (m q a)~ 2 log(m q a) since this would be a 
volume-dependent further contribution to the anomalous 
dimension of Q, and there can be no term in log(m 9 a) 
since the action is tree-level 0(a 2 ) improved [25]. A rig- 
orous proof of Eqn. (21 1 along the lines of Ref. [27] would, 
of course, be welcome. 

In the chiral limit m q — > 0, the term Wi that appears 

AQ) 

l 2,0 ■ 



on the right-hand side of Eqn. rt8| is a. 



E. Twisted spectral quantities 

The simplest spectral quantity that can be chosen 
within the framework of the twisted boundary conditions 
outlined above is the (renormalised) mass of the A me- 
son. The one-loop correction to the A meson mass is 
given by [T2] 



rn 



(i) 



-Zo(k) 



2m 



(0) 



(22) 



fc=(im™ ,0,m,0) 



where Zo(k) = 1 + [(ma) 4 ) is the residue of the pole of 
the tree-level gluon propagator at spatial momentum k, 
and is defined so that the momentum k is on-shell. 



Gauge invariance implies [B] 



a. 



a 

( m A,l) 

o 



(«»a,1) 
2,-2 

(m q a) 



0, 
. 



(23) 



The O (a s (ma) 2 ) contribution from improvement of the 
action is given by [12] 



>(!) 



(i) 



c 2 )(ma) z 



0((ma) 4 ) , (24) 



leading to the improvement condition 



A m A X) 

i 2,0 



(25) 



The next simplest independent spectral quantity is the 
scattering amplitude for A mesons at B meson thresh- 
old, which can be described by an effective AAB meson 
coupling constant A [13]: 

A = g ^Z(k)Z(p)Z(q) e, T x ^(k, V ,q) (26) 



where a twist factor of -^Tr([Tfc, r p ]r g ) has been factored 
out from from both sides, and the momenta and polari- 
sations of the incoming particles are (with r > defined 
such that E(q) = 0) 



k = (iE(k),k), 

p= H^(p),p), 

q= (0,q), 
e= (0,1,-1,0) 



k = (0, m, ir) 
p = (m, 0, ir) 
q = (— m, —m, — 2ir) 



(27) 



We expand Eqn. (26 1 perturbatively to one-loop order 



and find (up to 0(Jma) 4 ) corrections) 
1 



AW 
m 



1 



24 
1 - 



m 



fc dfc 7rn (fcj 



12 J dq 2 V g 



(?) 



fc =i-E(k) 

(28) 

90=0 



where I 1 ' 1 ) is the one-particle irreducible three-point 
function at one loop. The derivatives of the Feynman di- 
agrams contributing to the self-energy are computed an- 
alytically using automatic differentiation [2S] [30] . Con- 
tinuum calculations of the anomalous dimension and in- 
frared divergence give 



A\i) 

J Q.O 

x 2,-2 



1207T 2y 



The improvement contribution to A is |12] 

A imp — = 4(9c{ 1} - 7c { 2 1} )(ma) 2 + O ((ma) 4 
leading to the improvement condition 

4 (9c« - 7c«) = -4^ ■ 



(29) 



(30) 



(31) 



III. RESULTS 

To extract the improvement coefficients from our dia- 
grammatic calculations, we compute the diagrams for a 
number of different values of both L and m q with Nf = 1 , 
N = 3. At each value of m q , we then perform a fit in 
ma of the form given in Eqn. (19) to extract the coef- 
ficients a$ ,1 \m q a) for n 



a, 



(mA,l) 




(m q a) 



0, 2. Our fits confirm that 
0; an example is shown in Fig, jl] 



Performing a fit of the form in Eqns. ( 20|21 1 respec- 
tively on these coefficients, we are able to extract the 
analytically-known coefficients with high accuracy along 
with the required (ma) 2 contributions, as shown in Fig.pl 
Our results for the fermionic contributions are 



4™ A,1) = 0.00942(3) 
,(A,i) _ _o_3 52 (2) . 



a. 



2.(1 



(32) 



6 



A 



(0,*) 



V 



FIG. 3: The one-loop counterterm contribution to the static 
potential. 



Equating th ese res ults with the Wi of Eqn. we can 
solve Eqns. ( 25|31 1 for to obtain 



-0.025218(4) + 0.0110(3)^/ 
-0.004418(4) + 0.0016(3)^/ 



=S> 4 1} = 0.237088(46) - 0.1008(34)^/ (33) 

where the quenched (Nf — 0) results are taken from 
Ref. [T2] and we have propagated the errors by quadra- 



ture into c. 



(i) 



A. Effect on gauge action couplings 

The MILC and UKQCD Collaborations use a "tadpole 
improved" version of Eqn. (Ill, dividing each gauge link 
by a factor it . In addition, a factor of cq/uq is subsumed 
into the gauge coupling /3o = 6co/(3 2 Mq) that multiplies 
the plaquette term Pq 15J. The couplings multiplying 
the "planar rectangles" Pi and "parallelograms" P2 are 

mm 



fa 



(3o 
~20u 2 



( 



12tt 



.(i) 



V 5 



2«W In.. 



o 12^/3 (1) 
P2 = - 2 c\ 'a s 

5u 



(34) 



(with factors of 4-7T coming from converting from g 2 to 
a s ). The quenched radiative contributions have been an- 
alyzed in |15j and so we may write 



01 = 



ih 



0o 
20 

0o 



20m 2 , 



1 + 0.4805a s 



12vr 



(i) 



5 



■ 487rc^j. I n , 



0.033a s 



12tt 



-c (1) a 
2,/ ' 



(35) 



where now all the one-loop coefficients c^l contain only 
quark loop contributions. 

Plugging in the numbers for the HISQ action obtained 
in this work we find 



01 = 

02 = 



3o 



0.4805a s - 0.899(52)AT /Q : S 



20u 2 ^ 



~ [0.033a s - 0.0121(23)^/a s ] 

u Q 



(36) 



The full coefficient Co has here been absorbed into the 
gauge coupling, so the coefficient multiplying the pla- 
quette P is simply fa — fa. For the HISQ action, the 

fermionic contribution to cq (i.e. c^j in Eqn. ( 33 ) ) is 
large and sizeable shifts will be needed in fa to maintain 
a constant g 2 (or lattice spacing) as Nf is changed from 
(quenched) to Nf = 3 or 4. Whilst this is not a prob- 
lem in itself, it does make it more difficult to intuitively 
relate values of 0q to the lattice spacing. 

A more sensible choice is to absorb just the tree-level 
portion c °' = | into the gauge coupling. The overall 
gauge coupling is simply (3' Q — 10/(<7 2 Uq) . Using primes to 
denote couplings in this scheme, the coupling multiplying 
the plaquette in the action is now 



0' P 



0o 
0o [1 



Aire, 



(i) 



„(0) 



1.7876a s - 0.760(26)7V/a s 



(37) 



The remaining couplings in this scheme are 



0^_ 
20u 2 

P'o 
20u 2 

_0o_ 
20ul 



1 - 



(48^ Ci 1} 



2u 



(i) 



l + 2.2681a s 



(48.cS) 



[1 + 2.2681^- 1.659(46)^/^] , 



02 



l_2nfa (1) 
5u 2 ° 2 Us 

0'o 



[0.033a s - 0.0121 (23)^ 



(38) 



The factors multiplying the gauge coupling in fa an d fa 2 
are the same as this term is already 0(a s ). 



IV. RADIATIVE IMPROVEMENT AND THE 
STATIC POTENTIAL 

In this section, we seek to understand what effect the 
omission of the asqtad 0(Nfa s a 2 ) corrections to gauge 
action will have on physical observables measured in ex- 
isting non-perturbative Monte Carlo lattice simulations. 

As discussed above, the then-unknown 0(Nfa s a 2 ) 
contributions to Cj were omitted from the current gen- 
eration of three-flavour dynamical asqtad simulations. 
Perturbatively, this omission will lead to an imperfect 
cancellation of discretisation effects and a residual break- 
ing of rotational symmetry in the static quark potential. 
Similar effects are expected to be seen in the numeri- 
cal simulation results, and hence in the determinations 
of the scale-setting parameters ?i, fo derived from the 
static potential. 

Here we use one-loop perturbation theory to calculate 
the effect of the missing 0(Nfa s a 2 ) terms on the static 
potential and particularly on the scale-setting parameters 
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FIG. 4: The perturbative correction (top panel) and the 
correction for lack of rotational invariance (bottom panel) . 



fi, fn. The rationale for this and alternative approaches 
are discussed in Sees. Ill and fvl 



A. The correction to the static potential 



Including the missing 0(Nfa s a 2 ) corrections to the 
gauge action would change the lattice static-quark po- 
tential V(r) measured by the MILC Collaboration by an 
amount 5V(r). Here we estimate this change using one- 
loop perturbation theory. 

We do this by computing the spatial Fourier transform 
(for zero temporal momentum) of Go^(fc)A M „(fc)G„n(fc)j 
where G^ v (k) is the tree-level Symanzik improved gluon 
propagator at momentum k = (0, k), as shown in Fig. [3] 
is the 0(g 2 ) insertion into the gluon propa- 
gator arising from the perturbative expansion of the 
0(Nfa s a 2 ) corrections to the gauge action in Eqn. (11) 



with appropriate asqtad c^j couplings [BJ. Again, Feyn- 



man rules are derived using the HiPPy package [19ll20j . 

In all cases, the spatial Fourier transforms are carried 
out for a finite, periodic lattice of spatial volume L 3 . 

The gauge action has C3 = 0, but the static quark 
counterterm proportional to C4 is also omitted. This will 
also affect the success of the radiative improvement, but 
we do not consider its effect in this paper. 

The result for L = 20 is shown in the upper panel of 
Fig. |4j We expect the corrected lattice potential to be 
rotationally symmetric at 0(a s a 2 ), so the lack of rota- 
tional invariance in the upper panel of Fig. |4]is indicative 
of an equal and opposite breaking of rotational symme- 
try in the potential measured on ensembles that omit the 
(D(Nfa s a 2 ) radiative corrections. 



B. The effect of the correction 

To set the scale, the MILC Collaboration measure the 
static potential for a variety of on- and off-axis spatial 
separations: {V(ri)} with associated statistical errors 
{cr(rj)}. A least squares fit is performed using the fit 
function [3TJ &\ [33] 



Kfit(r) 
V coat (r) 

Vcorr(r) 



Vcottt(r) + 6sKorr(r) = E b ifj( r ) ' 

On h 2 r , 

r 

(V free (r) - i r < r cut 
1 otherwise 



which defines basis functions fj(r) with a vector of fit 
parameters b. V corr aims to account for the lack of rota- 
tional invariance at small r = |r|, with Vf ree the finite- 
sized lattice estimate for 1 fr from the Fourier transform 
of the (free) Symanzik gluon Coulomb propagator. We 
show this for L = 20 in the lower panel of Fig. [4j 

In more detail, the fit minimises the least-squared func- 
tion 



(39) 



We define the (weighted) average of operator A(r) over 
a set of measured r, as 




(40) 



The result of the least squares fitting is a vector of best-fit 
parameters b that obeys the linear equation 



where 



M jk b k = Xj b = M- 1 X (41) 



Mjk = (fjfk) , Xj = (fjV) • (42) 



Having done this, the lattice scale is set from the ana- 
lytic derivative of the V con t(r) function: 



: dV cont (r) 



dr 



, r ", • (43) 



We use "hats" here to stress that the scale parameters 
are measured in dimcnsionlcss lattice units. 

Two scales are commonly used: f\ from C\ = 1 (phys- 
ical value r*i = 0.317 fm [33]) and fo from Co = 1.65 
(physical value rn = 0.462 fm [33]). In general f\ is pre- 
ferred as the statistical errors on the static potential are 
smaller at shorter distances. 
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The corrected fits 



V. DISCUSSION 



Having calculated the 0(Nfa s a 2 ) corrections to the 
static potential, we can now calculate the effect of in- 
cluding SV on the best fit parameters bj, assuming that 
we carry out exactly the same fitting procedure as before. 

We would now minimise 

r' fh \ - \^ ( y ( r ') + a sSV(r t ) - V^jvi)) 2 ,,,, 



Given that b minimises L(b), we assert that &' = b + a s Sb 
minimises L'(b), with 

Sb = M~ 1 SX (45) 

and SXj = (fj SV). 

After finding 8b we can deduce the associated change in 
f\. We can either define this as Sf\ = fi(b + a s Sb) —fi(b) 
or, using a Taylor expansion of Eqn. d43|, 



(46) 



Sh _ If Sbi 5b 2 
The two methods give almost identical results. 

D. Results 



We looked at a range of ensembles listed in Table [T] 
using published values of uq to infer the strong coupling 
constant in the same way as MILC [33]: 



(47) 



4 log u 

Q, = . 

3.0684 

To estimate the effect the fermionic corrections would 
have on the scale setting parameters as measured by the 
MILC Collaboration, we adopt the same fitting function 
and we use the same fit range V5 < r < 7 and r cut = 3 
for the "fine" lattices and \[2 < r < 6 with r cut = 2.5 for 
the "coarse" and "very coarse" ensembles. 

We infer b\ and b 2 from published values for f\ and fo 
on given ensembles [32l [35] : 



1.65 



1 



b 2 r\ 



(48) 



For instance, on the (3 = 6.76, m u /m s = 0.01/0.05 coarse 
ensemble h = 2.60 (Ref. (32]), f = 3.76 (Ref. [35]) giv- 
ing 61 = 0.406, b 2 = 0.088. We then find 



Sri/h = -0.65%, Sfo/fo = -0.11%. 



(49) 



The shift in fi is larger because fi is smaller than f 
and SV is short-ranged. On the fine lattices, f 1 in lattice 
units is comparable to fo on the coarse lattices. The shift 
is therefore small. Results for other ensembles are given 
in Table [I] 

We have looked at various scenarios, e.g. different 
choices for the fitted range of {r^} and constraining some 
fit parameters to zero. Whilst the precise shifts do vary, 
the scale (and sign) of the shifts remain stable under such 
variations. 



Radiatively improved gluon actions are used in lat- 
tice simulations to give greater control over discretisa- 
tion effects and to reduce the uncertainty in continuum- 
extrapolated quantities. A typical example is the use of 
the Liischer-Weisz action in improved staggered simula- 
tions by the MILC and UKQCD Collaborations. 

We note that current unquenched simulations employ- 
ing lattice quark formulations other than improved stag- 
gered, such as domain wall or improved Wilson clover, 
generally do not use a radiatively improved action for 
the gluons; hence a calculation of the effects of fermion 
loops on the gluonic action is currently neither neces- 
sary nor useful for those simulations, but could readily 
be performed if and when simulations using such quark 
actions together with the Liischer-Weisz action will be 
undertaken. 

Simulations employing staggered quarks rely on the va- 
lidity of the "fourth root trick" , which has not yet been 
rigorously established. The purpose of this paper is not 
to engage in the debate about the validity of this pro- 
cedure, but merely point out that simulations using im- 
proved staggered quarks have produced results in excel- 
lent agreement with experiment so far. While we cannot 
completely discard the possibility that the observed scal- 
ing violation in <&b might be an indication of some more 
fundamental problem, we believe that our explanation 
for this scaling violation is more likely in the light of ex- 
isting evidence. In particular, we are able to replicate 
both the sign and the rough magnitude of the observed 
effect by a perturbative calculation. 

The shift in the radiative corrections due to the HISQ 
fermions in Eqns. ( 33|36 1 is surprisingly large, even com- 
pared to the coefficients for asqtad fermions [BJ. At 
first sight, this may seem like a surprise, since HISQ is 
supposed to be the more highly-improved action. How- 
ever, HISQ is designed to suppress taste-changing inter- 
actions coming from low momentum quark/high momen- 
tum gluon couplings, but the gluonic improvement coeffi- 
cients come from high momentum quark/low momentum 
gluon couplings, for whose suppression the HISQ action 
is not tuned. 

One consequence is that if the coefficient Co is sub- 
sumed into the leading factor of /3q, we expect to see 
large Nf dependent shifts in the value of f3o at fixed g 2 , 
and we also give results for an alternate scheme where 
only the tree level part of Co is included in the overall 
gauge coupling. 

The radiative corrections in the Liischer-Weisz action 
used by MILC in the asqtad simulations, however, omit 
the contribution from dynamical sea quarks. This con- 
tribution has recently been calculated at one-loop for 
Nf (massless) flavours of asqtad improved staggered 
fermions. The results are dramatic, leading to sign re- 
versals in some of the radiative coefficients. 

It is therefore conceivable that the omission of the 
0(NfCt s a 2 ) corrections leads to increased scaling viola- 
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TABLE I: MILC simulation parameters and shifts in scale setting parameters induced by omission of fermionic radiative 
corrections to the gluonic action. Smoothed fi values are from Ref. [331 136] . fo values are then inferred from the ratios fo/fi 
given in Ref. [33]. We have estimated uo for the very coarse ensemble. Lattice spacings are quoted as approximate guides; 
precise values may be inferred from setting n = 0.317 fm. 



Label 


a/fm (approx) 


L 3 


x T 


Sea quark 


h 


ro 


uo 


5f\/f\ 


Sfo/ro 










masses mi/m s 








(in %) 


(in %) 


very coarse 


0.18 


16 3 


x 48 


0.082/0.082 


1.805 (10) 


2.622 (28) 


0.8585 


-1.11 


-0.40 


coarse 


0.12 


20 a 


x 64 


0.02/0.05 


2.650 (8) 


3.828 (15) 


0.8688 


-0.63 


-0.11 










0.01/0.05 


2.610 (12) 


3.774 (20) 


0.8677 


-0.65 


-0.11 






24 3 


x 64 


0.005/0.05 


2.632 (13) 


3.834 (25) 


0.8678 


-0.64 


-0.10 


fine 


0.09 


28 a 


x 96 


0.0124/0.031 


3.711 (13) 


5.398 (28) 


0.8788 


0.01 


0.00 










0.0062/0.031 


3.684 (12) 


5.384 (27) 


0.8782 


0.01 


0.00 



tions in results from dynamical simulations when com- 
pared to quenched data. 

To properly establish whether this is the case would 
require a new set of dynamical Monte Carlo simulations, 
which is well beyond the scope of this study. An alterna- 
tive is to attempt a reweighting of the existing ensembles 
using factors e~ 5s based on the 0(Nfa s a 2 ) countert- 
erms. Such calculations notoriously suffer from a very 
poor overlap between the importance samplings of the 
original and reweighted ensembles for even minor changes 
in the action. This leads to very large statistical errors 
which will obscure any sought-for effect. 

We have therefore used instead one-loop lattice per- 
turbation theory to estimate the effect of the 0(Nfa s a 2 ) 
corrections on the static potential and the shifts in the 
scale-setting parameters f% and fo arising from the omis- 
sion of the fermionic radiative corrections for typical 
values of the simulations with 2+1 dynamical asqtad 
flavours. 

On fine (a ~ 0.09 fm) lattices, the shifts in f\ and 
fo are negligible (less than 0.1%) and will be at least as 
small on superfine lattices with a ~ 0.06 fm. On coarse 
lattices (a ~ 0.12 fm), omission of the corrections leads 
to fi being 0.6% too large, with fo unaffected. On very 
coarse lattices (a ~ 0.18 fm), f\ is 1.1% too large and fo 
0.4% too large. 

Overall, then, the omission of the 0(Nfa s a 2 ) leads to 
an underestimate of the lattice spacing on coarser lat- 
tices as defined using f\. Whilst numerically small, this 
effect is comparable to the statistical errors on a num- 
ber of quantities and therefore would lead to a measur- 
able increase in the statistical uncertainty of continuum- 
extrapolated lattice QCD predictions. 

Higher loop and non-perturbative effects will almost 
certainly change the exact value of the shift in fi, but 
are unlikely to alter our main conclusion: that the effect 



is measurable. 

Putting aside the static potential, an alternative ap- 
proach to fixing the lattice spacing is to use the 2S — IS 
mass splitting of T states. We have seen that that the 
correction 6V(r) is negative, so including the 0(Nfa s a 2 ) 
radiative corrections would decrease slightly both (lat- 
tice) T masses. Because 5V(r) is short-ranged, the effect 
on the 15 state will be larger than on the 2S* state since 
the IS wavefunction is larger at small r. The lattice mass 
splitting and thus the derived value of a will increase, and 
thus $b will get slightly smaller in such physical units 
on coarse lattices. To reliably deduce this fact from the 
T mass gap, however, we need to include the effect of the 
contact term, Eqn. d3J), which we do not yet know. 
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